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Abstract 

We investigate SU(2) gauge theory in a constant chromomagnetic field in 
three dimensions both in the continuum and on the lattice. Using a varia- 
tional method to stabilize the unstable modes, we evaluate the vacuum energy 
density in the one-loop approximation. We compare our theoretical results 
with the outcomes of the numerical simulations. 
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Typeset using REVT^ 



I. INTRODUCTION 



It was pointed out by different authors several years ago that for four- dimensional 
non Abelian gauge theories without matter fields in the one-loop approximation states with 
a constant chromomagnetic field lie below the perturbative ground state. This discovery 
stirred up a lot of interest. It was soon realized, however, by Nielsen and Olesen that 
these states are not stable due to long-range modes, the Nielsen-Olesen unstable modes. 
It was been, then, conjectured [^,^ that the vacuum field configurations which differ from 
the classical external field only in the unstable mode sector could stabilize the Nielsen- 
Olesen unstable modes. The conclusion was reached that first the constant chromomagnetic 
field would form a lattice, which, however, could be unstable under quantum fluctuations. 
These quantum fluctuations might result in the formation of a quantum liquid, the so-called 
Copenhagen vacuum 0. 

On the other hand, the problem of the unstable modes has been reconsidered from a 
different point of view 0. In Ref. henceforth referred as I, it has been showed that, 
by using variational techniques on a class of approximately gauge-invariant Gaussian wave 
functionals, the stabilization of the Nielsen-Olesen modes contributes to the energy density 
with a negative classical term which cancels the classical magnetic energy. Moreover the 
stabilization of the Nielsen-Olesen modes induces a further background fleld which behaves 
non analytically in the coupling constant and screens almost completely the external chro- 
momagnetic fleld. As a consequence, even in the strong-fleld regime, which is the naive 
perturbative regime, one deals with the non perturbative regime. This means that the cal- 
culation of the energy density even in the one- loop approximation is non perturbative 0. In 
other words the calculation of the vacuum energy is truly non perturbative and it mandates 
a completely non perturbative approach. 

In a previous paper we investigate four dimensional U(l) and SU(2) lattice gauge 
theories in an external constant background fleld. In particular, for the non Abelian case 
we looked for the effects of the unstable modes on the vacuum energy density. However, the 
conclusion was reached that the Monte Carlo data did not display effects due to unstable 
modes owing to flnite size of the lattice. Indeed, in four dimensions the unstable modes are 
long range. Then in lattice simulations one must approximate the continuum and, at the 
same time, work with a lattice big enough to allow unstable modes. So the external flelds 
which can be put on a lattice with periodic boundary conditions have strength P,p!0|: 



a^gB ~ (1.1) 



where n is an integer, L is the lattice size, and a the lattice spacing. 



On general grounds []IT|, one expects that it exsists a critical fleld above which the system 
becomes paramagnetic. As a consequence, we need external magnetic flelds lower than the 
critical fleld B,,. Using the estimation O] B^ ~ T^, Tc/Al — 40 and 



A,<.(,)./(,).(5^)-%.p{-5^i?) (1.2) 

we obtain for the minimal fleld Bmin which can be put on a lattice: 
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Be ~ L2 (40)2-^ 



(1.3) 



For instance, for /3 = 2.5 one need L > 12 to have Bmin/ Be < 1- Therefore, we see 
that in four dimensions working with external magnetic fields lower than the critical fields 
requires sizeable lattices. We can avoid these problems by considering three-dimensional 
gauge theories. Indeed, the three-dimensional SU(2) gauge theory is supcrrenormalizable. 
Moreover, as we shall discuss later, the unstable modes are not long range. 

The purpose of this paper is twofold. First, in Sect. II through Sect. V we consider 
the SU(2) gauge theory in a constant chromomagnetic field in the fixed-time Schrodinger 
representation. In Sect. II we show that the quadratic part of the Hamiltonian displays 
unstable modes. In Sect. Ill we discuss the stabilization of the unstable modes by using the 
variational procedure of I. Section IV is devoted to the calculation of the stabilized vacuum 
energy density in the one-loop approximation. In Sect. V we analyze the Nielsen- Ninomiya 
Ansatz. 

Second, in Sect. VI we reanalyze the problem by employing the lattice formulation of 
the gauge theories and compare the numerical results with theoretical expectations. Finally, 
our conclusions are drawn in Sect. VII. 



II. SU(2) HAMILTONIAN IN A BACKGROUND FIELD 

In this section we consider the (2-1-1) dimensional pure SU(2) gauge theory in the tem- 
poral gauge. We shall follow closely the method of I. In the temporal gauge — the 
Hamiltonian is 

H=ljd'x {{Et{x)f + {B%x)f} , (2.1) 



where 



and 



B'^ix) = le,,F^^{x) , (2.2) 



i^^.(f) = d,A';{x) - djA^ix) + ge'^''^At{x)At{x) . (2.3) 

Note that in two spatial dimensions the chromomagnetic field B°-{x) is a (pseudo)scalar. 

In the fixed-time Schrodinger representation the chromoelectric field Ef{x) acts as func- 
tional derivative 

on the physical states which are functionals obeying the Gauss' law: 

diS'^'^ + ge^'^'Atix)] j^MA) = . (2.5) 
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The effects of an external background field is incorporated by writing 



(2.6) 



where A1{x) is the background field, and r]°-{x) the fluctuating field. We are interested in a 
constant abelian chromomagnetic field. Thus we set 

A^{x) = 6''^6,2XiB , X = (xi, X2) . (2.7) 

It is straightforward to rewrite the Hamiltonian in terms of the fluctuating fields. We get: 



H = V— + + i7(3) + ^(4) 
2 



where 



6^ 



X 



6rif{x)6rif{x) 



+ vt{x)Dt {A;xH{x) - Df{xU{x) 



(2.8) 



(2.9) 



^^'^ = ?e,fee,vfc,e-''- / d'x 7^'{x)7^l{x)Df{x)4,{x) 



(2.10) 



H 



(4) _ 9_ 



^je,,erk'e''"'e^'''' J <fx ri]{xmx)4{x)4{S) 



(2.11) 



In Equations ( p.8[ - p.ll]) Df'{x) is the covariant derivative with respect to the background 
field: 



Df{x) = diS^" + ge^^'A^ix) . 



acb AC/ 



D1j{A,x) is the operator 



Dl\A,x) = -5,,Dl%x)Df{x) + 2gBe,,e 



Safe 



The Gauss' constraint 



can be rewritten as 



Dt{x) + ge^^'rit^x 



acb a/ 



6r]^{x) 







(2.12) 



(2.13) 



(2.14) 



As is well known |T^, the Gauss' constraint ensures that the physical states are invariant 
against time-independent gauge transformations. It follows, then, that the physical states 
are not normalizable. This means that we must fix the residual gauge invariance. Following 
I, we impose the covariant Coulomb constraint 



Dfix)vt{x) = . 



(2.15) 



As discussed in I, the functional measure in the scalar product between two physical states 
gets modified by the Fadeev- Popov determinant associated to the gauge- fixing ( |2.15D . 



We are interested in the vacuum energy. The one-loop approximation of the vacuum energy 
corresponds to consider the quadratic piece of the Hamiltonian. From equation ( p.8|) we get 



Ho = ^ + H('K (2.16) 



In the same approximation the Gauss' constraint reduces to 

A"'(f)||^0 (2,17) 

Equation ( |2.17| ) means that the physical states are functionals of the fields transverse with 
respect to the operator D°-^{x). 

To diagonalize it suffices to solve the eigenvalue equations 

05(f) = A0^(x) (2.18) 

with the conditions 

Dfix)<P^ix)=0. (2.19) 

The method to solve Eqs. ( |2.18| ) and ( |2.19| ) has been discussed in detail in Appendix B of 
I. So we merely write down the final results. We find the following transverse eigenvectors: 



N = k, k-e = 0, X{N,a = l) = P . (2.20) 



r,{N, a = 2; f) = j - J 0+(iV; f ) + -L [ +, | 07(iV; x) 











N={p,n) \{N,a = 2) = gB{2n + l), n>l 



where 



\{N, a = 




/ 1 


ai 




V2 


—i 
[ 


0+(iV,a 



(2.21) 



(j)+{N,a = 2;x) = ^ -i I F+{p,n + l;x) + +i I F+{p,n-l;x) 



( 1 



72 



V 



(2.22) 



In Equation (|2.22| ) «i = ^Jna and = ^/n + 1 a; the constant a is fixed by the normaliza- 
tion condition 



(fx (f)t{N, a; x) (jy^ (N', a; x) = S^n' ■ 



(2.23) 



Moreover 



^gB U + ^ 



(2.24) 



Finally 



V2 



-I |0+(iV;f) + 




V2 



/ 1 

V 



N = p, X{p^a = 3) = -gB , 



(2.25) 



with 



/ 1 



0+(p;f) 



V2 



-i F^{p, n = 0;x) 



\ 



. {p;x) = 0+(p;f) 



(2.26) 



The presence of the eigenvector (f)'}{N,a = 3;x) with negative eigenvalue implies that the 
quadratic Hamiltonian is not positive definite, and whence unbounded from below. Ac- 
cording we can decompose the fiuctuating field into a stable fiuctuation and an unstable 
one 



where 



r7f(f) =r/:,(f)+C(f) 



c{N,amN,a;x) 

Af,o=l,2 



(2.27) 



X 



N 



Whereupon, the quadratic Hamiltonian reads 



(2.28) 



(2.29) 



with positive definite. From Equation ( p.29|) it follows that the vacuum functional can 
be factorized as 



and the vacuum energy can be written as 



(2.30) 
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E{B) = V^ + E, + E^. (2.31) 

The contribution to the vacuum energy due to the stable modes is readily obtained in the 
one loop approximation. 

It is straightforward to show that the stable mode vacuum functional is: 

>v,(r^/) = exp|-^/ nn^^){^){gJ);H^Mn{^)^ , (2.32) 

with 

(a)J(x,y)= E 2A^(iV;a)0^(iV,a;f)0;*(iV,«;f). (2.33) 

N,a=l,2 

Then we obtain 

Es = - [ \k\ + V^—y^JgB{2n+l). (2.34) 

2 J (2Tif ' ' 27r V ; V ; 



In obtaining Equation ( p. 341) we used 

gB 



fdp\F{p,n-x)\^ = (2.35) 
J Ztx 



On the other hand, we have 



so that H^"* is unbounded from below. In order to stabilize H^'^ we must take into account 
the quartic coupling which is of order of g^. This means that we must take care of the 
Gauss' constraint at least at the order g"^ . 

The main strategy of these calculation has been discussed in Section III of I. In the next 
section we shall summarize the main achievements of our analysis. 



III. STABILIZING THE UNSTABLE MODES 



According to the previous section, we would like to stabilize the unstable modes by taking 
into account the quartic coupling. To do this, we must construct the variational vacuum 
functional which satisfies the Gauss' constraint up to order of g'^. We assume for the vacuum 
functional the form in Eq. (|2.30|) . The presence in Eq. ( p. 36 ) of a negative mass squared 
term suggests to try with a shifted Gaussian for the functional Ziirj): 



Z,{r],) = exp 



oo 



P§Pt k(§)-n^(§)] (^n)^U§,t) K|(t)-n[(t) 



(3.1) 



where 
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(GO"' = E 2p(p)0- (P, « = 3; (p, a = 3; f) 



(3.2) 



5:6(p)C(P,« = 3;x) 



(3.3) 



The vacuum functional satisfies the Gauss' law in the lowest approximation Eq. ( p. 17 ) 
In order to impose the full constraint Eq. ( |2.14| ), we write 



g',{7]) = exp [-,ir])] g,{r]) 



(3.4) 



and fix iteratively the functional ro(?7). 

To proceed, we evaluate the expectation value of the Hamiltonian on the vacuum func- 
tional ( p.4|) . Fortunately this rather long calculation parallels closely the one done in the 
case of three spatial dimensions. The reader interested in further details may consult section 
IV of I. 

In the one-loop approximation, the stable mode contribution to the vacuum energy has 
been already evaluated, Eq. ( p.34|) . For the unstable modes we get in the same approximation 



Eu = \ I d?x d?y 6{x - y) D'^^{A- x) 



G 



-1 



ab 



{x,y)+u1{x)u){y) 



+ <fx {G^)'-{x,x) + 



(3.5) 
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+ 



abc^ab c 



(fx ^u''j{x)ul{x)u'^',{x)u(,{x) 



_,\bb' ^ ^ ^ , ^ 11 

j _^{x,x) ul{x)ul,{x) + two permutations | . 



(3.6) 



In order to minimize Eu with respect to p{p) and b{p), we observe that from Eqs. ( p73| 
( ^ ) and it follows: 



u-'jix) = 



uf{x) 



g^{x) 



(3.7) 



with 



/27r V vr / 



(3.8) 



Using Equations ( p. 21 ) and (|3.7]), we recast into: 
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^ P 



P{P) 



gB 
pip) 



gB / d'^x g^{x)g (a;) + — / (fx g^{x)g [x] 



(3.9) 



As discussed in I, the necessary condition so that our configuration contributes to the energy 
density is: 

(3.10) 



(3.11) 



Sx g^{x)g {x) <x.V = L"^ . 

Moreover it is easy to show that the minimum of the energy is attained for 

g^{x)g'{x) = K 



where i^T is a positive constant. The condition (p.ll|) is fulfilled with the choice: 



h{p) = V2nKe'^P (3.12) 
in the thermodynamic limit L oo |13|. From Equations ( |3.11| ) and ( |3.12| ) we obtain 



^ P 



pip) - 



pip) 



- VKgB + ^VK^ 



+g^K Jd\ \F+{p,n = 0;x: 



(3.13) 



p pip) 

The last equation tells us that we can assume p{p) independent on p. Thus, by taking into 



account that 14 



gB 

27r 



and 



Y,\F+{p,n = 0; 



2 ^ gB^ 
2tx 



(3.14) 



(3.15) 



we get finally 



Eu ^g_B_ 
V 47r 



P- 



gB 



p 



Varying with respect to K, we get: 



K 



2gBK ^2 

g gBK + —K 

2tx p 2 



gB gBl 



5(2 271 p 



(3.16) 



(3.17) 



Inserting into Eq. (|3.16|) and neglecting terms of the same order as the two loop contributions, 
we obtain the simple result: 



V 



gB 

2 in 



P + 



gB 
p 



Now, the minimization with respect to p is straightforward: 



Whence 



P=9B. 



Eu B' {gBf (}^ 

V 2 27r V Gtt Gtt , 



(3.18) 



(3.19) 



(3.20) 



In conclusion we get for the total vacuum energy in the one-loop approximation the remark- 
able result: 



E{B) 
V 



B^ E^_B^ {gBf + O I'll 
2 V 2 2tt I Gtt Gtt , 



(3.21) 



El 

V 



(fk 

(2^ 



n=l 



(3.22) 



The last two terms in Equation ( p. 21 ) are the contributions due to the stabilized unstable 
modes. Note that, as in three spatial dimensions, the unstable modes contribute to the 
vacuum energy density with a negative classical term which cancels out the classical magnetic 
energy term. 



IV. THE STABILIZED VACUA 



In the previous section we evaluated the vacuum energy density in the one-loop approx- 
imation. If we neglect the unstable mode contribution our result should coincide with the 
real part of the one-loop effective potential: 



(Pk 

W) 



gB 



(4.1) 



n=l 



Obviously V^'^\B) is divergent. However the theory is superrenormalizable. This means 
that once we subtract the free vacuum energy density we are left with a finite result. As we 
show below this is the case and the final result 



Ar(2)(5) = V^^\B) - \/(2)(o) = ^ - -^^^ 



1 



V2 



47r ^ V2 



(4.2) 



where C,{z) is the Riemann's Zeta function, is in agreement with the calculation by H. D. 
Trottier p[. 

Indeed, by using the identity: 
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°° ds d 



we have 



gB /■°° ds d 



v^V^ '^■s 



-g_B(2n+l)s 



27r JO y^y/s ds 



^+gBs 



-gB, 



Moreover, we observe that 



/ 



d'^k 



\k\ 



I ds d /I 



(27r)2 An Jo v^y^ ds \s 

Thus, Equations ( [4 .41) and (|4.5| ) ensure that 

3 r d^k 



/ 



F(2)(o) 



Whence 



^ ^ 2 27r 271 



1 



A change of the integration variable recasts Eq. ( [4 .71 ) into ||T6 

f"^ dx d 







2 


2n 


B^ 




~2 


271 



1 + 



1- v^C 



n^/x dx I 

1 r 

2' 2 



1 

2x 



where C{Z)(l) is the generalized Riemann's Zeta function. 

y/2-l 



Observing that ITB 



47r 



C 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



we obtain the promised result Eq. 
The one-loop effective potential ([4.8|) has a negative minumum which, however, lies in a 
region where the one-loop approximation is not trustworthy. 



When we take into account the unstable modes we get |]T7 
AE{B) , V2- 1 



V 



+ 



87r2 



■c 



3/2 



o 



(4.10) 



In Equation ( [4.10| ) the coefficient in front of {gB)^/'^ is positive. As a consequence, the 
minimum of the vacuum energy density is attained for 5 = 0. 
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It should be stressed that the calculation of the full contribution to the vacuum energy 
is non trivial even in the one- loop approximation. As a matter of fact, Equation ( [4.10 ) does 



not include the contributions to the energy arising from the interaction of the stable mode 
with the induced background field u'^ix) which behaves non analytically in the coupling 
constant. Nevertheless, the cancellation of the classical magnetic energy term due to the 
stabilization of the unstable modes is sound. So that we can write in general 

^ . , O (ill?) . ,4.11, 

In order to avoid the unphysical situation where the vacuum energy density decreases without 
bound for increasing external magnetic field, the constant 03/2 should be positive, so that 
the minimum is again B = 0. 

It is also interesting to note that the total background field is now given by: 

A'^ix) = A{x) +ir{x) . (4.12) 
From Equations (|3.7|), ( ^.111 ), and ( ^.17] ) we obtain for the full field strength tensor 



4 = ±^gB (5,i5,2 - 6,26,1) . (4.13) 

On the other hand the external field strength is 

F^f = d,A] - d,Al = B (5a5,2 - 5,26,1) . (4.14) 



It follows that 



P3 
-^12 

F'cxt 
12 




(4.15) 



Equation ( [4.15[ ) can be compared with the result we should have obtained by neglecting the 



unstable modes: 



|l = 1 ■ (4-16) 



Equation ( ^.15| ) tells us that the background field induced by the stabilization of the unstable 
modes strongly screens the external magnetic field. 

In conclusions, the main results of this section are summarized in Eqs. ( 4.12| ) and ( 4.151 ). 



Note, however, that our results Eqs. (|4.12|) and ( [4.15|) are valid only in the thermodynamic 



limit L — i> 00. In view of the comparison with Monte Carlo simulations on a finite lattice, 
it is worthwhile to discuss another Ansatz due to H. B. Nielsen and N. Ninomiya which 
can be realized even on a finite lattice. 
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V. THE NIELSEN AND NINOMIYA ANSATZ 



In this Section we discuss the Ansatz by Nielsen and Ninomiya 



9 [X) 



1 



1/4 



/27r V TT / 



(5.1) 



(5.2) 



with if TV constant to be fixed by minimizing the vacuum energy density. Using the definition 
of the third Jacobi 6 function |jT6[ 



^3(^,g) = l + 2Eg"'cos2 



nz 



n>l 



we rewrite Eq.( |5.1|) as 



where 



q = e 



(5.3) 



(5.4) 



(5.5) 



z = - {X2 + ixi^ 



It is easy to show that 

J (fx g+{x)g-{x) = Kn ' J dxidx2 e-^^"' \e^{z,q){' 

'gB^'l^ 



K 



N 



TT 



dxi dx2 



and 



d^x 



g [x)g [X) 
Inserting into Eq.(p.9|), we get 



2 g_B 

271 



^NiT' [^3(0, g)]^ / dxidx2 . 



p{v) 



gB_ 
p{p) 



gB 



2n 



2tx 



(5.6) 



(5.7) 



(5.8) 



9BK^ Viz)" y^^^^l i^-^] \m. V 



+ a'^Y^ J l^'^iP^ n = 0]x)\ g'^{x)g (x) 



(5.9) 
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Assuming p{p) independent on p, and using Eqs. ( p.l4| ) and ( p.l5|) , we recast Eq. (|5.9|) into: 



V An 



gB 

P 



gB 



gB\ 

271 J 



1/2 

) + 



g^gB 

2 271 



^ 2 9B fgB\'/' 1 



Varying with respect to Kj^i we obtain 



271 



1/2 



gB 



gB 



g'[9siO,qW 27r [^3(0,?)]^' 



whereupon 



Eu^g_B_ 

V A7T 



P + 



gBe 



B' 



with 



[^^3(0,g)]^ 



2[^3(0,g)]^ 



1 . 



Note that ^3(0, g) = 1 + 26"^ + 26"^^ + ■ • ■ ^ 1.0864, so that e > 0. 
The minimization with respect to p is now straightforward. We get: 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



P = ygBe , 



(5.14) 



- ^ r2 , V^r^D^3/2 



(5.15) 



Finally, the total vacuum energy density is 

1 



AE{B) _ 1 
~ 2 



V 



1 - 



+ O 



B'- 



(gB) 



3/2 



27T 



4^^y 



(5.16) 



We stress once again that the coefficient of (gBY^'^ in Eq. ( |5.16D does not include the 
contributions due to the interaction between stable modes and the induced background field 
<(x). 

Note that the coefficient of the {gBY^'^ term is negative, so that Eq. ( 5.16| ) displays 
a negative minimum. However, the minimum lies in a region where two-loop terms are 
sizeable. 

It is instructive to evaluate the field-strength tensor. The only non-zero component of F^^, 
is Ff2- We obtain 
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FUx)=B-gK^ i^^y\~a^^i\e,iz,q)\' , (5.17) 

one can show that the chromomagnetic field forms a square lattice in the xi — X2 plane 

with lattice constant ajsf = qB. 
Defining a space averaging 

Fl) = , (5.18) 



/ d?x 

and using Eq.( |5.7| ), we get 



Fl,) = B-gK^ f^V^' . (5.19) 



\2ti J 

Finally, Equations ( |5 . 1 1|) and (|5.14| ) bring on 

<l!l=il 1 ^ I 1 1 (5 20) 

Ff^ \ [^^3(0,g)Py'^27rV^[^3(0,g)]2 ^ 

where = B. 

The Nielsen and Ninomiya Ansatz is interesting in view of the comparison with the 
lattice approach to be discussed in the next section. Indeed the chromomagnetic lattice 
could be more easily realized on the space-time lattice if the two lattices are commensurate. 
In the lattice approach the external chromomagnetic field is quantized due to the periodic 
boundary conditions: 

a gB = —n , n integer (5.21) 

where a is the space-time lattice spacing, and L the linear lattice size in lattice units. From 
Equation (|5.21|) we get 

^ (5.22) 



a V2n 

Thus the lattices are commensurate if 2n is the square of an integer. 



VI. BACKGROUND FIELDS ON THE LATTICE 

From the previous section we learned that the calculation of the vacuum energy is truly 
non-perturbative even in the one-loop approximation. Thus, we need a non-perturbative 
approach which can be furnished by the lattice formulation of gauge theories. In Ref. |^ we 
studied the four- dimensional pure SU(2) lattice gauge theory in an external magnetic field. 
We looked for the effects due to unstable modes. Due to the limited size of the lattice we 
used in the Monte Carlo simulations, we failed in observing the unstable mode effects. The 
reason for such a failure have been already discussed in the Introduction. 
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The situation should improve considerably working with three-dimensional SU(2) gauge 
theory. Indeed in c? = 3 we can perform Monte Carlo simulations on lattices of considerable 
size. In addition, in two spatial dimensions the unstable modes are not long range. 
The pure SU(2) gauge theory is implemented on the lattice through the standard Wilson 
action 



X,fl>U 



X) 



(6.1) 



where U^^{x) is the elementary plaquette in the (/i, i/)— plane at the lattice site x. 
The external background field on the lattice can be introduced via an external current |]8|,p!0 
For the reader convenience, we briefly summarize our method. 
In the Euclidean continuum the background action reads 



S 



B 



d'x j:xx)a:{x) 



Using the classical field equations: 



r(x)=D:\A)F'{x 



-Fj^ being the field strength tensor built from the external fields A"(x), and 



S^'F^fix) 



we get from Eq. (|6] 



d,Al{x) 



d.A^ix) 



(6.2) 



(6.3) 



(6.4) 



(6.5) 



To discretize the action Eq. (|675|), we must define the Abelian-like piece d^A^{x) — d^A^^^x) 
from the plaquette variables. To do this, we use the so-called Abelian projection [p^|l9| . 



To implement the Abelian projection, first we fix the gauge by diagonalizing an operator 
X{x) which transforms according to the adjoint representation of the gauge group: 



V{x)X{x)V\x) = diagonal matrix. 
After that, we rewrite the gauge- fixed links 

U^{x) = V{x)U^{x)V\x + ^x) 



(6.6) 



(6.7) 



as 



where 



U,{x) = W,ix)U^ix) 



U^{x) = diag 



(6.8) 



(6.9) 



arg • 



{[u,ix)U . 



(6.10) 
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The Abelian-like links (|6.9| ) are called the Abelian projection of U^{x). The Abelian pro- 
jected plaquettes t/^j, are built form the Abelian projected links in the usual manner. Clearly 
we have 



diag 



It is, thus, natural to define the Abelian field strength tensor as 



(3 tr 



(3 sm9^,{x 



As a consequence, we are led to consider the following background action 



S 



B 



(6.11) 



(6.12) 



(6.13) 



Taking into account the periodic boundary conditions of the lattice, we write 



/3sinC(^) 



(6.14) 



2tt 

l2 



(6.15) 



the n^f^^{xys being integers. In Equation ( |6.15D , L is the lattice size. 

We are interested in a constant Abelian chromomagnetic field. In this case the lattice action 
reads: 



S = +[3 U,J 

X,fl>U 



x] 



/5^sin^^f sin^ 



12 



X 



(6.16) 



The quantity of interest is the vacuum energy density at zero temperature in presence of the 
external magnetic field. In the naive continuum limit one can show easily that the vacuum 
energy density is given by 



E [F^f 



with 



P.(Frr) = l-itrf/i2 (Fff) 



(6.17) 



(6.18) 



E [l-^trt/3.(Fff) 



1=1.2 



(6.19) 



We performed Monte Carlo simulations with the action (|6.16|) on lattices of size L = 20, 
and L = 40 in the weak coupling region (3 > 7. We measured the quantity 



AE (F^f ) = E (F^f ) - F(0) = 

= p{{Ps{Fu')) - {Psm + (F.(o)) - {Pt{F^r))} ■ 



(6.20) 
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Our simulations were done in the unfixed gauge, which corresponds to set V{x) = U in 
Eq. Q. 

We stress that the background action Eq. ( |6.13| ) depends on the gauge-fixing. On the 
other hand, in the continuum hmit the energy density is gauge invariant. This means that 
one should verify that the lattice definition ( |6.17| ) does not depend on the gauge-fixing 
procedure. In our previous study in four- dimensions we found a very weak dependence on 
the gauge-fixing. Moreover we are interested in the effects due to the unstable modes. As 
we have already discussed, the unstable modes modify in a dramatic way the vacuum energy 
density. So that the effects we are looking for should manifest even without any gauge fixing. 

After discarding 500 sweeps, we collect 500 measurements (1 every 5 sweeps). Statistical 
errors were evaluated by the jackknife algorithm. 

In order to reduce the statistical errors the difference {Pg — Pt) was evaluated directly during 
Monte Carlo runs. Note that the vacuum energy difference can be also measured directly 
by observing that, within statistical errors, we have 



tr U: 



3i 



cxt 
12 



(tr f/3.(0)) = (tr U^m 



(6.21) 



In Figure 1 we display the adimensional energy density AE {F^2^) j versus the adimen- 
sional external field strenght j for four different values of /?. 

In general the vacuum energy can be affected by finite lattice effects. To avoid lattice 
artefacts we used quite large lattices. We found that for 20^ lattices there are no sizeable 
finite size effects up to /5 = 10 (compare full and open circles in Fig. 1). On the other hand 
for /? > 10 we need larger lattices (see full and open triangles in Fig. 1). As a consequence 
for /5 = 12 and 15 we doubled the lattice size, so that we feel that our numerical results are 
reliable. 

In the continuum the adimensional energy density can depend on the unique adimensional 
combination at our hands, namely B j g^ . So we can write 



AE{B) 



f 



B 



5" \r/ 

In the one-loop approximation we have 



where 



0-3/2 



1 

2^ 



l-^C(3/2) 



if we neglect the unstable modes. Including the unstable modes we obtained 

V2-1 



«3/2 



■am 



02 = 



in the approximation of section |IV|, and 



«3/2 



1 
2^ 



a2 



[^3(0, g)]^ J 



(6.22) 



(6.23) 



(6.24) 



(6.25) 



(6.26) 
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for the Nielsen and Ninomiya Ansatz. Note that the unstable modes cause a drastic reduction 
of the classical magnetic energy term. 

On the lattice there is another dimensional quantity, namely the lattice spacing. Whereupon 
Equation (|6.22| ) is a non trivial check for the Monte Carlo outcomes. Any deviations from 



the scaling law ( |6.22| ) can be considered as an indication of the granularity of the lattice. 



A glance at Fig. 1 shows that the scaling law (|6.22|) is satisfied quite well in the weak field 
strenght region F^2^/g^ < 1. But, for strong field strenght F^2^/g^ > 1 we observe a small 
deviation from the scaling law (|6.22|) . This found a quite natural explanation. Indeed, in 



order to approximate the continuum, the magnetic length 1/ y^'-Fif* should be greater than 
the lattice spacing. Thus, at fixed lattice spacing, this condition gets worse by increasing 
the external field strength. 

In Figure 1 we compare our numerical results with the theoretical expectations. The 
dotted line is the one-loop effective potential (|6.24| ), the dashed line is Eq. ( p.25| ), and the 
dot-dashed line Eq. (|6.26|) . Note that the theoretical calculations are restricted to the one- 
loop approximation. However this is not a limitation at all. Indeed the two-loop terms are 
important for 

g^gB > jgBf/' 
27r 271 ^ 27r ' 

i.e. 

B 



< 



g3 47]-2 



(6.27) 



Unfortunately our Monte Carlo simulations do not allow us to appreciate the vacuum energy 
differences already for B/g^ ^ 0.4 . 

From Figure 1 it is evident that the Monte Carlo data strongly disagree with the one- 
loop effective potential Eq ( |6.24D . Indeed the vacuum energy density is about an order 



of magnitude smaller than the classical magnetic energy. We feel that this reduction can 
be ascribed to the unstable modes. As a matter of fact, our numerical results agree with 
Eq. ( |6.25| ) for F^f' / g^ < 1 . On the other hand for F^f' / g^ > 1 the data can fitted by using 



Eq. ( |6.23|) with 02 ~ 10 ^. However, as we have already stressed, the complete cancellation 



of the classical magnetic energy can be attained for configurations which satisfy Eq. ( p.lll ). 



This condition can be fulfilled only in the thermodynamic limit. When one deals with a 
finite lattice, it is natural to expect the presence of a small classical-like term in the vacuum 
energy density. 

In order to have a further check, we looked at the Abelian and non Abelian chromo- 
magnetic field strength. Indeed, a further clear signature of the unstable modes resides in a 
drastic reduction of the non Abelian chromomagnetic field with respect to the Abelian one. 
In Figure 2 we display the ratio (f/^) / Fff versus FfgVs'^- From Figure 2 we see that 
the system feels a rather strong Abelian chromomagnetic field. Note, however, that the 
Abelian chromomagnetic field is strongly affected by lattice artifacts. In Ref. ||2T]] we found 
a similar situation even in the U(l) lattice gauge theory. The situation is quite different for 
non Abelian chromomagnetic field. As Figure 3 shows, the ratio (-^^2) /-^i2* depends only 
on Fl2^/g^ (the small deviations in the strong field strength region have been already dis- 
cussed). Moreover the non Abelian chromomagnetic field strength is reduced by one order 
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of magnitude with respect to the Abehan one. In Figure 3 we also compare the numerical 
data with Eqs. (|4.15|) and ( ^.20 ). A few comments are in order. Without unstable modes 



and in the one-loop approximation we should have obtained 



Fext 
1' 



1 . (6.28) 



12 



Equation (|6.28|) is the dotted line in Fig. 3. We see that Eq. (|6.28|) is in total disagreement 



with the numerical results. In addition, if we take into account that a small quadratic term 
in the energy density adds a constant term of the same order in Eq. (|4.15| ), then we can 
conclude that the data corroborate Eq. ( |4.15| ) at least in the region Fff / > 0.5. As 



concerns the weak field strength region, on the one hand we expect that in Eq. (^4.15|) there 
are sizeable effects due to higher order contributions. On the other hand, the numerical 
results are not reliable because the magnetic length became comparable to the lattice size. 
Thus we are confident that Figs. 1 and 3 together give evidence of the unstable modes on 
the lattice. 

Let us conclude this section by comparing our approach to the one recently proposed by 
H. D. Trottier and R. M. Voloshyn These authors consider three-dimensional lattice 
gauge theories in a background field. The background field is induced by means of an 
external current Eq.( |6.2| ). Thus in the continuum the background action they used coincides 
with Eq. (|6.5| ). For the U(l) gauge theory the action of Ref. p2[ agrees with our lattice 



background action Eq. (|6.13|) . However the authors of Ref. p2| do not enforce periodic 
boundary conditions, so that the external background field strength is not restricted by the 
constraint Eq. (|6.15 ) [^]. As concerns the SU(2) gauge theory, Trottier and Woloshyn adopt 



a different discretization of the Abelian field strength tensor: 

FU^) = (U,, - [U,, U,]]^ . (6.29) 

It should be stressed that, whatever discretization one chooses the total action is no longer 
SU(2)-invariant, but only U(l)-invariant. This means there is always the freedom of reducing 
the original local SU(2) invarance to a local U(l) invariance with a suitable gauge fixing. In 
particular, the authors of Ref. used the discretization Eq. (|6.29|) without gauge-fixing. 



As concerns the numerical results, the authors of Ref. found that the vacuum energy 
density is more than an order of magnitude smaller than the classical magnetic energy over 
the whole range of the applied external magnetic field, in qualitative agreement with our 
results. 

However, for ^ 1 the vacuum energy density of Ref. is negative. This result 

looks puzzling to us because we do not expect that in the weak external field strength region 
the vacuum energy density displays a dramatic dependence on the discretization of the 
Abelian field strenght tensor. To check this point, we performed Monte Carlo simulations 
for F^2^/g^ < 1 by using Eq. (|6.29|) in the background action. It turns out that the vacuum 
energy density difference is positive even though it is about a factor two smaller than our 
previous results. This disagreement can be ascribed to the different discretization adopted 
for the Abelian field strength. It turns out, however, that the discretization Eq. ( |6.29| ) leads 
to a vacuum energy density which is close to Eq. ( |6.25| ). Thus the method of Ref. p2 



corroborate the evidence of unstable modes on the lattice. 
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Contrarily to our findings the authors of Ref. p2[ found that the vacuum energy density 
is negative for Fff^/^f^ < 1 . We do not beheve that this disagreement can be ascribed to the 
different boundary conditions. We do not yet understand the reasons of this discrepancy. 
We hope to clarify this point in future studies. 



VII. CONCLUSIONS 

Let us conclude by summarizing the main achievements of this paper. We investigated 
the three-dimensional SU(2) gauge theory in a constant chromomagnetic field, both in the 
continuum and on the lattice. As in the four-dimensional case, we found the presence 
of unstable modes. The stabilization procedure of the unstable modes introduces a new 
background field with non analytic behaviour in the coupling constant. The most striking 
consequences of the induced background field are the cancellation (or a drastic reduction) 
of the classical magnetic energy, and an almost complete screening of the applied magnetic 
field. Moreover, we found that the non analytic behaviour of the induced background field 
makes the calculation of the ground-state energy highly non trivial. We studied, then, the 
problem with the non perturbative techniques offered us by the lattice approach to the gauge 
theories. We feel that we have convincingly put out the evidence of the unstable modes. 

We would like to end by briefly discussing some consequences of our finding. As is known 
since long time, in the non Abelian gauge theories in the perturbative one-loop approxima- 
tion, states with a constant chromomagnetic field lie below the perturbative ground state. 
It follows that these states could be a better approximation to the true ground state than 
the perturbative one. In particular, it was believed that from these states one could 
set up an approximate ground state which confines color charges. However, the presence 
of unstable modes leads to a muddled state. In the case of (2+l)-dimensions, we argued 
(see the discussion after Eq. ( |4.10| )) that, after the stabihzation of the unstable modes, the 
states with a constant chromomagnetic background field are not energetically favoured with 
respect to the the perturbative ground state. This means that these states are not relevant 
to the confinement at least in three spacee-time dimensions. However, if we belive that the 
confinement mechanism does not depend on the space-time dimensions, then we are led to 
the same conclusions even in the more interesting four-dimensional case. 
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namely the vacuum energy density was close to the classical value | (-P'lf*)^- 
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FIGURES 



FIG. 1. Vacuum energy density versus the applied magnetic field. Full and open symbols refer 
to L=20 and L=40 respectively. Squares correspond to /? = 7, circles to (3 = 10, triangles to 
f3 = 12, and diamonds to /? = 15. The curves are discussed in the text. 



FIG. 2. The Abelian chromomagnetic field strength versus the applied magnetic field (symbols 
as in Fig. 1). 



FIG. 3. The full chromomagnetic field strength versus the applied magnetic field (symbols as 
in Fig. 1). The dashed line is Eq. ( 4.15| ), the dot-dashed line Eq. ( 5.20| ), and the dotted line is 
Eq. f|6:2F 
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